Incommensurate Charge Order Phase in Fe20BC>3 due to Geometrical Frustration 
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The temperature dependence of charge order in Fe20B03 was investigated by resistivity and 
differential scanning calorimetry measurements, Mossbauer spectroscopy, and synchrotron X-ray 
scattering, revealing an intermediate phase between room temperature and 340 K, characterized by 
coexisting mobile and immobile carriers, and by incommensurate superstructure modulations with 
temperature- dependent propagation vector (|,0,r). The incommensurate modulations arise from 
specific anti-phase boundaries with low energy cost due to geometrical charge frustration. 

PACS numbers: 71.30.+h, 64.70.Rh, 61.10.Eq, 76.70.+y 
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Geometrical frustration, the inability of interactions 
to be simultaneously satisfied due to the geometry of the 
underlying crystal lattice, in the case of magnetic inter- 
actions is known to be responsible for incommensurate 
spin modulations [1] and for various exotic ground states 
[2[. The Coulomb-interactions leading to charge order 
(CO) Si can also be geometrically frustrated [5] . How- 
ever, the effects of geometrical charge frustration are less 
well understood, despite of the classical CO example Q 
occurring on a frustrated lattice. Observations of incom- 
mensurate charge modulations are typically attributed 
to charge-density- wave-like Fermi surface nesting scenar- 
ios [6|, [3], but unusually complex CO patterns observed 
in spinels "devil's staircase" -type CO modulations 
found in NaV^Os under pressure [9j], and the appearance 
of a new type of ferroelectricity upon charge ordering in 
LuFe204 [10], may be related to charge frustration. 

The warwickite crystal structure of Fe2 0B03 (Fig. [T] 
inset) is prone to frustration, and we recently observed 
a CO superstructure with integer valence states (clearly 
distinct from a classical charge-density-wave) at low T 
[Tlj . Here, we present a detailed study of the thermal 
evolution of the CO in Fe20B03, establishing the pres- 
ence of an unanticipated intermediate phase [12[ with 
coexisting mobile and immobile carriers. Surprisingly, X- 
ray scattering indicates that the superstructure is incom- 
mensurate in the intermediate phase with a T dependent 
propagation vector (~,0, r). We show how the incom- 
mensuration arises from geometrical charge frustration 
via the proliferation of low-energy anti-phase-boundaries, 
as proposed also for chemical order in binary alloys [l3[ . 
Similar incommensurate phases may occur in other frus- 
trated systems with "binary order" (e.g. Ising spins) and 
thus our results not only provide the first example of in- 
commensurate phases with ionic (and therefore binary) 
CO, but have implications for the broader study of geo- 



metrical frustration. In particular, we propose Fe2 0B03 
as an ideal system to study ordering dynamics on frus- 
trated lattices. 
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FIG. 1: (Color online) Bulk physical properties of Fe2 0B03 
indicating two phase transitions, a— d: T dependence of resis- 
tivity p along a (a), heat flow of differential scanning calorime- 
try (b), estimated monoclinic angle (3 (c), and key parameters 
describing the Mossbauer spectra (d): fraction of the contri- 
bution to the spectra by Fe ions with no discernible electron 
hopping (■, left) and hopping frequency v of the remainder 
(O) right). Inset in a: Crystal structure at 355 K (after [111]). 
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We studied single crystals from the same batch as those 
studied in [ll|. The electrical resistivity was measured 
with standard four-probe geometry, using gold-sputtered 
contacts. Differential scanning calorimetry (DSC) was 
performed, with a Perkin-Elmer Diamond DSC, on sev- 
eral powdered crystals, with temperature sweep-rates be- 
tween 10 and 50 K/min giving consistent results. The de- 
tails of Mossbauer spectroscopy and electronic structure 
(GGA+U) calculations are as in [Tlj . 

Synchrotron X-ray scattering was performed on two 
crystals with a mosaic spread of 0.20(8)° at the MU- 
CAT sector 6IDB of the Advanced Photon Source in the 
Argonne National Laboratory, with the photon energy 
tuned to 7.05 keV. All scattered intensities shown are 
normalized by an ion chamber monitor to account for 
changes in the incoming beam intensity. Due to the crys- 
tal mosaic spread and the presence of twinning, we deter- 
mined the monoclinic distortion by rotating the crystal 
to the nominal (1,3,3) position and scanning the detec- 
tor angle (similar to a powder diffraction experiment). 
Using the a, 6, and c values obtained in [11] the mon- 
oclinic angle /3 was estimated, with error bars given by 
the uncertainty in a, 6, and c. A further study was per- 
formed at 6IDD with photon energy ~ 98keV and an 
image plate system, using a rocking technique as in [14J . 

Between 250 and 400 K, resistivity p(T) and DSC data 
(Fig. QJi,b) show two separate, well defined phase transi- 
tions on cooling (/warming) at 340 and 280(/308)K, the 
former of which corresponds (Fig. []J) to the monoclinic- 
orthorhombic structural transition. The transitions de- 
lineate low, intermediate, and high T phases. Ac specific 
heat had not revealed the lower transition [ll| , which we 
attribute to its insensitivity to irreversible contributions 
with hysteresis larger than the ac amplitude. In contrast, 
DSC is sensitive to all contributions to the specific heat. 
The estimated total entropy change through both tran- 
sitions (~0.6 J/mol/K/Fe ion) still is only about 10% of 
the configurational entropy associated with perfect CO. 
Two likely reasons for this will be discussed below. 

Mossbauer spectra (Fig. [2j) indicate divalent and triva- 
lent Fe ions distributed over two structural sites, as be- 
fore [15[. However, whereas there is no discernible hop- 
ping of electrons between Fe ions at low T (hopping be- 
low ~ 1 MHz is not resolvable), the description of the 
high T spectra requires hopping with a T dependent fre- 
quency v (Fig.[T]i, o), which remains quite low even above 
400 K [16|] . In the intermediate T phase between the two 
transitions, e.g. at 325 K, contributions with and without 
electron hopping coexist (Fig. [TJi, ■). Isomer shifts and 
quadrupole splittings are equal for both contributions, 
and follow a standard T dependence. 

We further characterized the three phases by X-ray 
scattering using synchrotron radiation. In the low T 
phase, weak additional reflections (Fig. [3fci) indexing to 
(h+7},k,t) are observed, as expected [11]. Superstruc- 
ture reflections are sharp along h and k. Along t they 
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FIG. 2: (Color online) Mossbauer spectra (+) at selected T on 
powdered Fe20B03 crystals, together with fits (black lines). 
Different contributions to the fit are indicated by colored lines. 
Bright blue and magenta: contributions from Fe sites with no 
electron hopping (two doublets for each color correspond to 
two distinct Fe valences) . Blue and red: contributions from Fe 
sites with non-zero electron hopping. Red and magenta dou- 
blets, with the larger quadrupole splitting, may be assigned 
to the sites on the outer chains I, IV and the blue and bright 
blue doublets to the sites on the inner chains II,III |15( . 



are >25x broader than fundamental reflections, and in- 
dicate a correlation length, obtained by the inverse of 
the full width at half maximum, of only £ c ~ 12 c. This 
suggests differently ordered domains (likely up and down 
diagonal CO [11], see Fig. (4^l) of size large along a and 
6, but small along c. These micro-domains correspond 
to an imperfect overall CO, which is one reason for the 
lower-than-expected entropy. In the intermediate phase 
superstructure reflections are sharp along t and are split, 
corresponding to an incommensurate modulation, with 
propagation vector (^, 0, r), and r(T) varying from about 
0.4 at 340 K to 0.2 at 280 K (Fig. [3]o). Above 340 K a 
long-range ordered superstructure no longer exists, but 
very weak and broad reflections with (h + \, £-\- \) in- 
dex (Fig. [3]i) indicate persistent short-range correlated 
fluctuations, with correlations mainly along a, the chain- 
direction. The correlations are not static, but the dynam- 
ics is relatively slow (Fig.[T]i). This dynamic short-range- 
order likely also contributes to the "missing entropy" . To 
test the incommensurability of the modulations in the in- 
termediate temperature phase, data were collected in 1 K 
intervals with high-energy X-rays. No indications for any 
"lock-in" to commensurate values are visible in Fig. [3b, 
thus the modulations are truly incommensurate rather 
than forming a "devil's staircase" as in [9[. 
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FIG. 3: (Color online) Commensurate and incommensurate 
superstructure reflections from synchrotron X-ray scattering. 
Fundamental reflections at integer (h,k,£) typically have peak 
intensities in excess of 10 4 counts/monitor, a: Scattered in- 
tensity at (|,fc,£) and 100 K, showing one of the superstruc- 
ture peaks corresponding to the doubled unit cell [1J| . Black 
line: half- maximum- intensity contour, b: Intensity at (§, 3, i) 
vs T. An estimated beam heating of 5K was corrected. Inner 
white ticks indicate T at which data were taken, c: Intensity 
after background subtraction at (|,0,£) measured every IK 
in a different experimental setup with an image plate system, 
d: h, k, and £ scans through (f ,3, f) at 350 K. 



To elucidate the micro- domain formation at low T, and 
the incommensurate CO at intermediate T, we consider 
the energy of various CO configurations. The electro- 
static energy is minimized by having as few same- valence 
nearest neighbors as possible. This requires strictly al- 
ternating valences along a within each chain and estab- 
lishes "ordered chains" as the basic CO unit, consistent 
with the correlations in the fluctuation regime (Fig. [3jl). 
Geometrically, the 2D lattice of chains consists of two 
sublattices offset from each other by ~ (1.6 A): One sub- 
lattice comprises chains numbered I and III in Fig. [4] and 
Fig. [T] inset, the other chains II and IV. The electrostatic 
interactions between any two chains not belonging to the 
same sublattice are geometrically frustrated (Fig. EK). 
Therefore, flipping all valences in all chains (changing 
their phase) of one sublattice does not change the energy. 
Thus, the two configurations of Fig. [4^ are degenerate in 
energy, as are the two of Fig. [4)3. This twofold degener- 
acy is lifted by a monoclinic distortion (ribbon tilting), 
and hence the system will distort, thereby gaining energy 
However, according to our GGA+U calculations the 
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FIG. 4: (Color online) Different charge order configurations 
within the four-chain ribbons of the Fe2 0B03 structure (Fig. 
□ inset). Fe 3+ and Fe 2+ are drawn as filled and open circles, 
respectively, a: Diagonal charge order, the ground state con- 
figuration. Degeneracy of the two configurations leads to do- 
main formation, with opposite sense of monoclinic distortion 
(llj . b: Alternative "zig-zag" CO with only slightly higher en- 
ergy than the diagonal groundstate. All other configurations 
are intermediate between diagonal and zig-zag. 
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FIG. 5: (Color online) Impact of geometrical frustration on 
the charge order. a,b: Phase-relationship between two charge- 
ordered chains: For chains offset by f (a), inter-chain electro- 
static interactions are frustrated, and two configurations are 
degenerate in energy. Without offset (b), the anti-phase con- 
figuration has lower energy, c: be projection of the Fe2 0B03 
structure. The arrangement of chains leads to geometrical 
frustration even for interactions (thick lines) between chains 
within one sublattice, with possible phase relationships as in 
b (the other chains are crossed out). Dashed line: possible 
anti-phase-boundary (see text). 



energy gain associated with this distortion is very small, 
-0.2% of the overall CO gain for diagonal CO (Fig. UK), 
insufficient to prevent the facile formation of domains. 

Neglecting the small monoclinic distortion there is no 
net interaction between the two sublattices so that fur- 
ther considerations can be restricted to one of them. Any 
two chains of the sublattice are either in-phase or anti- 
phase (Fig. [5t>) with a lower electrostatic energy if they 
are anti-phase. Nevertheless, this does not remove all 
frustration because there is additional geometrical frus- 
tration even within one sublattice, due to the arrange- 
ment of chains projected in the be plane (Fig. [5b): several 
"frustration triangles" are readily visible (highlighted in 
red, green, blue). Although this frustration is not exact 
(triangle sides have unequal length), various configura- 
tions are expected to be very close in energy. Indeed, 
while our GGA+U calculations established a diagonal 
CO (Fig.Hk) ground state they also show that, e.g., 
the total energy gain due to zigzag CO (Fig. 0t>) is only 
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1.9% [18] smaller than the one due to diagonal CO, de- 
spite having more chains in-phase rather than anti-phase. 

Due to this near-degeneracy of configurations with dif- 
ferent phase relationships, resulting from a hierarchy of 
geometrical frustration, the energy cost of creating an 
anti-phase boundary in the ab plane (a possible location 
is indicated in Fig. [5t) between up and down diagonal 
CO domains (Fig. [4^) is small, explaining the small £ c . 
In contrast, there is no geometrical frustration within a 
chain, leading to a much larger £ a . £5 is large also due 
to the association of different CO domains with mono- 
clinic distortions of opposite sense. A boundary in the 
ac plane between domains of different monoclinic angle (3 
is not possible (short of cleaving the crystal), because the 
lattices will not match. In contrast, there is no problem 
matching domains with different (3 in ab or be planes. 

Incommensurate CO in classical charge-density- wave 
systems [i| and colossal magnetoresistance manganites 
[3] has been associated with Fermi-surface nesting. How- 
ever, because Fe20B03 is far from metallic (Fig. Hk,d), 
a nesting scenario can not account for the incommen- 
surate modulations in the intermediate phase. A bet- 
ter analogy to our case of essentially perfect ionic CO 
is found in the chemical order of ions in binary alloys. 
In fec-structured alloys, such as Cu-Au, incommensu- 
rate phases between chemically ordered (commensurate) 
and disordered phases are often observed, and can be ex- 
plained as a modulation with statistical period arising 
from anti-phase boundaries easily created due to degen- 
eracy and frustration of the bonds on triangular networks 
fl3| | . Given that in the intermediate phase of Fe20B03 
the incommensurate modulation is in the c direction and 
the anti-phase boundaries _L c are the ones with a small 
energy cost, we deduce that a similar scenario applies. 

The strong and continuous T dependence of the prop- 
agation vector then implies that the boundaries move 
easily in the intermediate phase, and thermal excitations 
should cause them to fluctuate. Fluctuating boundaries 
reduce the monoclinic distortion (as observed, Fig. [Tfc) 
and they imply that a fraction of Fe ions, those near the 
boundaries, change their valences. Fe ions change their 
valence by electron hopping, providing a natural expla- 
nation for the presence, in the intermediate T phase, of 
both contributions with and without electron hopping 
in the Mossbauer spectra (Fig. [T]i). The coexistence of 
mobile and immobile electrons, reminiscent of nano-scale 
inhomogeneities prominent in manganites and cuprates 
[Hj], but coherently ordered in Fe20B03, further stabi- 
lizes the intermediate phase as a compromise between CO 
and delocalizing tendencies influenced by the geometri- 
cal frustration (without which the transition to a charge 
ordered state likely would occur at much higher T [20[). 
Similar compromise phases may occur in other correlated 
electron systems with no chemical disorder. 

In summary, we discovered an intermediate T phase 
of CO in Fe2 0B03, which is characterized by the coex- 



istence of mobile and immobile carriers, and by an in- 
commensurate superstructure, the latter of which can be 
explained by specific anti-phase boundaries created easily 
due to geometrical charge frustration, similar to incom- 
mensurate chemical order . It should be interesting to fur- 
ther scrutinize the dynamics |5[ of this model-system of 
incommensuration from geometrical charge frustration. 
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